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Introduction 

In quantum field theory, the principle of locality is enforced by an exact light cone. Whenever two (bosonic) 
observables are spacelike separated they have to commute, so that neither can have any causal influence on the 
other. In nonrelativistic quantum mechanics, no explicit request for locality is imposed and it is possible in principle 
to signal between arbitrarily distant points in an arbitrarily short time, as can be seen from simple perturbative 
arguments. Nevertheless, the simple perturbation analysis shows that such an influence must decay exponentially 
with the distance between the observables. The seminal work by Lieb and Robinson 1 1 1 has made the preceding 
statement rigorous. In essence, it states that any quantum system whose Hilbert space is composed of a tensor 
product of local Hilbert spaces and whose Hamiltonian is the sum of local operators will have an effective maximum 
speed of signals. The lightcone is only approximate since it possesses an exponentially decaying tail. 

Recently, the Lieb-Robinson bounds have received renewed interest and have been applied to many areas in the- 
oretical condensed matter and quantum information theory. In particular, they were used to prove that gapped 
systems have exponential clustering in their ground state [2, 3J. The structure of correlations in quantum lattice 
systems is of particular importance in condensed-matter theory and quantum information theory [4-8 , 26|. In |7|, it 
is shown that the creation of entanglement can proceed with a speed that is limited by the Lieb-Robinson speed and 
it is shown that this is useful for proving general scaling laws for entanglement. In [5] these techniques have been 
exploited to characterise the creation of topological order. Other developments of significant interest include ll9l [T0] , 
where it is shown that it is possible to entangle macroscopically separated nano-electromechanical oscillators in an 
oscillator chain and that the resulting entanglement is robust to decoherence. To do this, the coupling between the 
oscillators is first turned off and the system is brought to its ground state. Once in the ground state, the coupling 
between the oscillators is turned back on non-adiabatically so as to produce a squeezed state which propagates en- 
tanglement over macroscopic distances at a speed related to the Lieb-Robinson speed. Such a system would be of 
great interest for its possible application as a quantum channel and as a tool for investigating the boundary between 
the classical and the quantum world. 

A systematic application of this technique can yield results for quantum spin models in two dimensions whose 
exact solution leads to intractable mathematical complexities. In a domain where the vast majority of results are of 
computational nature, analytical results are very important. 

One area of condensed matter physics in which Lieb-Robinson bounds have yet to leave their mark is quantum 
phase transitions. Quantum phase transitions i TTHT5] are phase transitions which occur in the quantum regime 
when the interaction strengths are varied (as opposed to thermodynamic quantities like temperature and pressure 
in the classical case). So far, with the exception of [16, 17] (in which only two interactions have been considered), the 
Lieb-Robinson bounds have been independent of the relative interaction strengths. When it comes to quantum phase 
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transitions however, the relative interaction strengths are of the essence. In future work, we would be interested in 
applying the techniques developed in this paper to Quantum Graphity[18 19], a model of quantum gravity where it 
is believed that space emerges from spacetime via a quantum phase transition. 

In this paper, we start by reviewing the Lieb-Robinson bound in Section [I] As a precursor, we show in Section 
III] how the dependence on the interaction strength can be obtained for a system with two interactions. This was 
done earlier in [16], but the derivation is presented here in greater detail. This is also done in lTl7l for the cases of 
unbounded operators. We then apply this formalism to the Ising model which is a system of two interactions. In 
Sectionlllllwe go on to generalise the previous results to the case of a finite number of interactions and apply it to the 
XY-modeL 



I. THE LIEB-ROBINSON BOUND 



In its simplest form, the Lieb-Robinson bound can be derived in the following way. We suppose that the total 
Hilbert space is the tensor product of local Hilbert spaces: 

H = ®H X . (1) 

X 

Each local Hilbert space % x is associated to an edge or vertex of a graph G with vertices having a maximum valency 
v. The total graph is thus associated to the total Hilbert space. On this Hilbert space we have a local Hamiltonian 
which defines the evolution: 

Hl oca l = ( 2 ) 
XcG 

where the operators $x are R-local. We define an operator to be R-local if it has support on a Hilbert space with 
an associated subgraph X of G such that the diameter of X in graph distance is less than R and such that X n Y = 
=$> $y] = 0. The norm of these local operators is also assumed to be less than some positive number K. The 
Lieb-Robinson bound is a bound on the commutator of two local observables separated by a distance d at a time t 
such as 

f(t):=[O P (t),O Q (0)}, (3) 

where the operators Op and Oq have support on the Hilbert spaces associated with the subgraphs P and Q which 
are separated by a graph distance d. Taking the derivative with respect to t and rearranging using the Jacobi identity, 
we have that 

f(t) = -i ([/(t),^(t)] + [Op(i),[$x(t),O Q (0)]]) (4) 
xeZp 

where Zp contains subgraphs of diameter at most R which have a non-empty intersection with P. Taking the norm 
and integrating we obtain 

\\[Op(t),O Q (0)]\\<\\[Op,O Q }\\+2\\O P \\ /* ds\\[$ x (s),O Q (0)}\\. (5) 

The preceding equation can then be iterated to give 

\\[Op(t),O Q (0)}\\ < \\Op\\\\O q \\ J2 M^a„ (6) 



n=0 

where 

n 

a « : = E • • E IIII^^q" <K n l n e^ nR -^ (7) 
x 2 ez Xl x n+1 £Z Xn j=i 

for some positive number A and where we have set X x .= P. We define 8* = lifXnF^0 and Sy = otherwise. 
We also define / := maxx |Z(X)| < oo for diameter(X) < R. A bound on a n is obtained by bounding ||3>Xj II by K 
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and noticing that there are at most I terms in each sum such that Sq — and nR < d. Using this bound for a„ we 
obtain a form of the Lieb-Robinson bound 



\\[O P (t),O Q (0)}\\ < \\0 P \\\\0 Q \\exp(-X(d-v\t\)), 



(8) 



where v < (2Ke XR )/X. The preceding inequality implies that the "speed of signals" in the system is effectively less 
than v. Anything travelling faster than this speed will be exponentially suppressed. 

From the preceding derivation of the Lieb-Robinson bound, two limitations immediately surface. First, passing 
from Eq.||4]| to Eq.|[5]|, the norm of the local operators is used; hence it is necessary that the local operators composing 
the Hamiltonian be bounded. Second, it is evident from using ||$Xj|| < K in Eq.d7l that in obtaining the Lieb- 
Robinson bound, one has ignored the interplay between different types of operators (different types of interactions) 
weakening the bound in the case where not all of the local operators have norm K. In other words, the bound can 
be tightened by classifying the different operators making up the Hamiltonian into sets of "interactions" where all 
operators in one set commute with each other and come with the same coupling constant in the Hamiltonian and 
then keeping track of how they appear in the bound. Though the the first limitation cannot be removed entirely 1 29 1, 
it was partially removed in |17J for a class of unbounded operators having bounded commutators. The purpose of 
this paper is to alleviate the second limitation. 

Before we proceed, however, we point out a stunning implication of the existence of the Lieb-Robinson bound: A 
relativistic light-cone structure emerges in the continuum limit of certain local (non-relativistic) quantum systems. 
More precisely, the continuum limit of any quantum system, for which the Lieb-Robinson bound is applicable, will 
exhibit a sharp emergent light-cone structure with a finite signaling speed. A heuristic argument can be given by 
replacing graph distances with metric distances in some units. For instance, let A be the unit graph distance in some 
standard unit such as metres or centimetres. If we now rewrite the Lieb-Robinson bound using this standard unit 
instead of the graph distance, we get 



where the Lieb-Robinson speed, vlr, is measured in the standard units of length per unit time and the distance 
d(P, Q) is measured in the same units of distance. It is now clear that if we take the continuum limit of the theory, 
A — > 0, Eq.Q will be exactly zero if d(P, Q) > vlr \t\. A complete analysis of the emergence of a strict light cone in 
the continuum limit has been shown in [20 1 for general harmonic quantum systems. 

The Lieb-Robinson bound allows us to conclude that any discrete R-local quantum theory with a Hamiltonian 
consisting of uniformly bounded or commutator-bounded operators will have an exact emerging light-cone structure 
in the continuum limit. Since many theories of quantum gravity are built from discrete structures such as graphs, the 
Lieb-Robinson bound could provide a natural mechanism for the emergence of relativity in those theories without 
requiring local Lorentz invariance a priori. 



Before taking into account the subtle interplay between a finite number of interaction types in the Hamiltonian, 
we will derive the Lieb-Robinson bound for a Hamiltonian with two interactions. This gives us the intuition needed 
to derive the bound for a system with more than two interactions. We consider R-Local quantum systems with 
Hamiltonians composed of bounded operators containing only two interactions. We will begin by defining precisely 
what an R-local quantum system is. 



R-local quantum systems are quantum systems with interactions of finite range R. More rigorously, let G be a 
graph (finite or infinite) having a maximum valence such that to every edge e we associate a the Hilbert space H e 




(9) 



II. THEORIES OF TWO INTERACTIONS 



A. R-Local Quantum Systems 



1 This ensures that the continuum limit of G be locally compact. The requirement will be needed so that the sum 
encounter later, converges for all real w. 



^-r , which we will 
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and to every vertex v we associate a Hilbert space H v . The total Hilbert space of the system is then defined to be 

Htot = (g) n e ® (g) ft„. (10) 

e an edge of G v a vertex of G 

We then take the evolution of the system to be i?-local for R, a natural number. By this we mean that the evolution 
is governed by a Hamiltonian composed of a sum of operators, each having support on a Hilbert space which is 
associated to a subgraph of G having a diameter less than R (i.e., the graph distance between any two points of the 
subgraph for which the operator acts is at most R). By association, we also refer to the diameter of an operator as the 
diameter of its associated subgraph. Also, we use the notation that if $ is an operator that does not have an explicit 
time dependence, then it is considered to be an operator evaluated at time t = 0, i.e., $ = $(0); furthermore, we 
have = e itH $e~ itH . 



B. General Hamiltonians With Two Bounded Interactions 



A common type of Hamiltonian is one with two coupling constants and two different types of operators. We say 
a Hamiltonian is a general Hamiltonian with two coupling constants if it is of the following type: 

ies jeSi 

with the properties that for a e Z 2 = {0, 1} and G S% we have [<&*, = and ||$* || = \\$° a \\ = 1. In other 
words the operators are normalised and operators of the same type at different locations on the graph commute. 
Here, the set S a is the set of labels of the subgraphs of G for which the operators of type a have support. 
Let T(q, m) be the subgraph associated with the operator <i>™ and let (a, b) 6 Z|. We define the function 

in analogy with Eq. (|3}. We will generalise the Lieb-Robinson bound as calculated in lTl6l to R-local Hamiltonians 
of bounded operators having two interactions by deriving a bound for K l a \(t). Note that since we are dealing with 
bounded operators, it is possible to absorb the norm of the operators into the coupling constants and we may thus 
assume \\^\\ — 1. In the following, we will show that if P and Q are two separate regions of the graph which are 
small, as compared to the typical size of the total graph, and Op and Oq are bounded local operators acting on the 
subgraphs P and Q respectively, then 

\\[O P (t),O Q (0)]\\ <(7exp^(2| ev / M^i-d(P,Q))), (13) 

where C is a positive constant depending on Op, Oq, and H. Here 7 and £ are positive constants depending on the 
graph structure, and d(P, Q) is the graph distance between P and Q. This substantiates the claim |5J that information 
cannot propagate at speeds greater than Vlr = 2^ey7iofti without suffering exponential suppression of the band- 
width as a function of distance. We will call this limit on the speed of information propagation the Lieb-Robinson 
speed. We obtain this result through recursion, as described in what follows. 



Taking the derivative of Eq.l 12 1 with respect to t, we get 



= -ih a+1 [[&®,*%i(t)]M] (14) 
-ih a+1 J2 Kw>*a^iW] + Kw,[^iW-^ 



iK+i E ^iW)] + HWi) E KW.K+iW-^]]. (is) 



where the Jacobi identity was used to obtain the third equality and for i € S q with q € Z 2 , we define Zi to be the set 



of k € S q +i such that T(q, i) n T(q + 1, k) 7^ 0. Note that if a e Z 2 and a = 1 then a + 1 = 0. Thus integrating Eq.l 15 1 
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we get 



K i a \ j (t)=T 2 (t)K i a \ j (0)-ih a+1 fds J2 Rt^KiW-^J], (16) 



where T 2 (t) is a unitary evolution given by the unitary matrix U 2 (t) = exp ( — ih a+ i J2 i2<£Z . ®a + i{t)) with T 2 (t)0 = 
U2(ty OU 2 (t) for any operator O. Taking the norm of Eq.l 16 1 we obtain 

iixlv'(*)ii<ii^V(o)ii+2ft. +1 f ds e (it) 

= ||K(*) ) ^]ll+2^ + i||^ 1 (0)ll fds ]T \\K\i(s),H}l 
<2Sf+2h a+1 fds ]T ||[< 2 +1 (s),^]||, 



2Sy+2h a+1 fds Y, WK+i b (s)l (18) 



where we define 



S k = f i if r(o<,*) n r(a fc , *) ^ 0, (19) 

5 1 otherwise, 

and where the last inequality is obtained by realizing that || <f>|J | < 2 and is necessarily zero if T(c, A:) and T(d, I) 
have an empty intersection. Applying Eq.{l8} to itself, we get by induction that 



K b j (t)\\<Mj: {2ltlV ^ r ( E - E 4")- (20) 



where M = max j \J~^, \f^\- Note that c„ counts the number of chains associated to n local operators (and thus n 
subgraphs T(c, k)) starting with T(a, ii) and ending with T(b, j). In all such chains a T(0, k) is followed by a I) 
where T(0, k) n T(l, an d vice versa. If we repeat exactly the same procedure, but this time with Op(t) instead 

of ^ (t) and Oq(0) instead of where Op, Oq are operators defined on some regions P, Q we get 



|[Op(t),Og(0)]||<iV||O F ||||O Q ||Mxf; ^ 2 -^g a ^ ( £ ... £ (21) 

n=0 " i2eZ ;i i„ + ieZi, 



where A/ is defined as previously described and iV is the (finite) number of local operators of the Hamiltonian 
intersecting P. Here c n counts the number of chains associated to n local operators (and thus n subgraphs T(c, k)), 
starting with T(a, ii), which must intersect P, and ending with T(b, i n ), which must intersect Q. Furthermore, we 
can always find a bound of the following type for c„: 

Cn <M7 n e A( f- d) , (22) 

where A is an arbitrary positive real number. This is because the T(a, i)'s have a diameter of R or less. Hence, if the 
distance d between the initial and final point is greater than nR then there are no possible chains of n local operators 
linking the initial and final point. Since there are always fewer than v other local operators intersecting T(a, i) (y 
must be finite because all vertices have a valence of not more than Q, there is a maximum of v n possible chains of n 
local operators starting from any given position. Thus, we certainly have that 

c n < v n e x{Rn - d) . (23) 
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Using Eq.d22l in Eq.|21 1 we obtain 



\\[O P (t),O Q (0)]\\ <J&ex P A(2 V / Mi'^eti-d), (24) 

where M — M MN\\ Op \\ \\Oq ||. This gives us an upper bound on the speed information can travel 

v LR ^2 v ^h 1 ^, (25) 

which can be minimised by choosing A = £. We thus get the following upper bound on the speed information can 
travel 

v LR = 2-eyfhohl. (26) 

C. Lieb-Robinson Bound on a Lattice 



In the case where the graph is a regular lattice, it is possible to give a bound of the form given in Eq.(22l as a 
function of the lattice and the Hamiltonian. We show here how this can be done. 

For the Hamiltonian given in Eq.jTT) on a regular lattice we have the following properties (in addition to the 
properties on a general graph): By passing from a $o to a <£>i in an operator chain, one always moves along the lattice 
by a lattice distance of D\ and by passing from a $i to a <&o in an operator chain, one invariably moves along the 
lattice by a lattice distance of Dq; in addition, for a given i, the finite number n a ^ of j's such that T(a, i) and T(b, j) 
have a non-empty intersection is independent of i and j and depends only on a and b, that is, the interaction types. 
Since the chains must alternate between operators of type and operators of type 1, the number of operator chains 
of length n starting at a specified point on the lattice, but with an unconstrained end point, is, at most, 

n -n In^o 



c„<max , ,./ } ^n ^ini^o . (27) 

I V "-i^o V n o^i J 

where the max prefactor takes care of both even and odd chains. In this case, for d(P, Q) > (D +Di ), there cannot 

be any chains of n operators linking the initial and final point. This means that for A > 0, e x {^i~( Do+Dl )~ d ( p 'Q}) is 
necessarily greater than 1 if there exists at least one operator chain of length n linking P and Q. Thus, instead of 
Eq.([27| we can write 

c n < MV^^~o n e X ^ (Do+Dl) - d{P ' Q) \ (28) 
where M = max j J*j=^, \f^^\ e #( D o+- D i), Referring back to Eq.|22l we obtain: 



7 = A/no-^ini^o, (29) 

( ^ <30) 

Combining these results with Eq.([26) gives 

vlr = evAio^ini^o^oMA) + -Di)- (31) 

Of course, the restriction to lattices is a significant one; however, as we see in what follows, it is possible to obtain 
an equivalent bound on c„ for more general graphs. 



D. Lieb-Robinson Bound on a Homogeneous and Isotropic Graph 



Let us now suppose that we have a general graph, and not necessarily a lattice. In this case, the D's of Section II C 
gain a dependence on the specific local operators $' and that is, instead of simply Do and D\, we have D 1 ^ 3 and 
D{^ k . Similarly, the n's gain a dependence on i: Instead of simply rip^i and ni^o, we have n Q _ > _ 1 and n{^ Q . For a 
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given operator chain of n operators , 4% +1 , • ■ ■ , 4> b n 1 , <frb+i}' tne en d points of the chain can be separated by, at 

most, a graph distance of YJCl D^ ih+1 . Hence, if d(P, Q) > YZ=l D^ ikJtl , the chain cannot link P to Q. 

Let us consider chains of n operators starting at i\. Let %k and jk be different chains such that i\ = j\. Then if we 
have 



n-l 



1 



chain i k=i 



n-l 

V 

n ^ 
chainj fe=i 



Oft 



' Jfc+1 and, 



\ 



n 



chain i fc=i 



n-l 

\ chain j fc=i 



(32) 



(33) 



and so the Lieb-Robinson speed will be the same in all directions. By w we mean that the quantities converge for 
large n and are thus independent of the path the chain makes on the lattice. Moreover, if the quantities in Eqs. 
{32} and | [33| depend on i n/ then the Lieb-Robinson speed will vary in different regions of the graph. Finally, if the 
quantities in Eqs. |32| and |33| depend on n, then the speed will depend on the distance d. 

We say a graph isnomogeneous if the quantities in Eqs. p2| and < [33) do not depend on n. Furthermore, we say a 
graph is isotropic at i\ if Eqs. (32 1 and d33l do not depend on i n , in which case the Lieb-Robinson speed will be the 
same in all directions. If these conditions are satisfied, the Lieb-Robinson speed will be constant in all directions and 
areas of the graph (or subgraph) considered. 

It follows that in the case of a homogeneous and isotropic graph, we can define 



\ 



n 



D = 



1 



chain i k=i 

n-l 

Enifc-s-'ifc+i 

chain i k=i 



(34) 



(35) 



for some chain of length n where n is large. By definition, these quantities will not depend on the particular chain 
chosen. We then have that n n is equal to the number of chains of n operators starting at a given point and ending 
anywhere. Note that for d(P, Q) > nD there cannot be any chains of n operators linking the initial and final points. 

This means that for A > 0, e A ( nD ~ d ( p > < 2)) is necessarily greater than 1 if there exists at least one operator chain of 
length n linking P and Q. Thus, we can write the analogue of Eq.|28| by 



(^\(nD-d(P,Q)) 



c n < n n e 
which gives 

7 = n, 

Combining these results with Eq.|[26) gives us 
v L r = 2eDn\/h hi. 



(36) 

(37) 
(38) 

(39) 



As in Sections II C we require a bound for the coefficient c n in order to provide a Lieb-Robinson speed in terms of 
the parameters of the system. Let 




chain ifc=i o fc -^o fc +i) 



- jipo + Di), 



D 



n 2^chain i fc=i c ' k 



lattice, 
h&i graph, 

lattice, 
h&i graph, 



(40) 



where "h&i" stands for homogeneous and isotropic. Using the definitions of Eq.|40l, we have that Eq.d39J is the 



general equation for the Lieb-Robinson speed for both lattice systems and homogeneous and isotropic systems. 
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E. Example: The Ising Model 

The quantum Ising model describes a spin network with a preferred alignment that is subjected to a magnetic 
field transverse to this preferred direction. In one dimension, the model is exactly solvable and is known to exhibit 
a second-order phase transition at zero temperature 1 14J . The quantum Ising chain has a wide range of applications 
from condensed-matter physics [21 ] to quantum gravity ||22| . 

The quantum Ising chain is an example of a bounded system on a lattice having a Hamiltonian with two coupling 
constants. Accordingly, we can use the framework developed previously to calculate the Lieb-Robinson bound. 

The Hamiltonian for the quantum Ising chain is given by 

i 

where a x and a z are Pauli matrices. At zero temperature the Ising chain exhibits two phases: paramagnetic (g > 1) 
and ferromagnetic (g < 1) with a quantum critical point at g = 1. By varying g, that is, the transverse magnetic field, 
the system can undergo a quantum phase transition from the paramagnetic phase (disordered) to the ferromagnetic 
phase (ordered) or vice versa. 

In order to calculate the Lieb-Robinson speed we need an upper bound on the coefficients c n defined in Eq.|[22|. 
The c„ coefficients represent the number of paths of length n consisting of alternating erf and af a", 1 operators such 
that the subgraphs corresponding to the operators intersect (in this case only nearest neighbours). Since we are 



dealing with a lattice we can use Eq.|29l and Eq.p Ot. There are 2 possible steps from a vertex to an edge and 2 



possible steps from an edge to a vertexTThus by Eq.|29 1, in this case, 7 = v2 x 2 = 2. Furthermore, by passing from 
a vertex to an edge, and vice versa, we move by half ofa lattice spacing, thus by virtue of Eq.|[30), £ = 2. 

It follows from Eq. p6) that the Lieb-Robinson bound for the speed of propagation in the quantum Ising model is 
found to be 



vlr = 2ey/gJ. (42) 

Close to the quantum phase transition (which happens at g = 1), at small deviations from criticality, \g — 1| -C 1, 
it is possible to take the continuum limit which gives a relativistic Majorana fermion field theory in one dimension 
with a speed of light of c = J and a mass of m = (g — 1) J [? ]. This means at up to a constant of order 1, 2e, which 
we expect due to the approximations made, the Lieb-Robinson gives the right value for the speed of information 
propagation, at least near the phase transition where it can be checked. 



III. THEORIES OF FINITELY MANY INTERACTIONS 



Now that we have established in detail how to calculate the Lieb-Robinson bound for two interactions, we will 
apply the intuition we gained to derive a bound for any finite number of interactions. The essential intuition which 
we gained for two interactions was that, in an operator chain, because one type of interaction commutes with inter- 
actions of the same type, each operator in a chain must be of a different type than the preceding one. In the case of 
only two types of interaction, this requirement uniquely defines the interaction type of each operator in the chain; in 
other words, they must alternate. However, in the case of rn interactions, the situation is more complicated because 
we have a choice of as many as rn — 1 interaction types for each operator along the chain. The way we will deal with 
this added complication is by considering periodic chains of operators. Since the bound for an arbitrary chain will 
be the geometric mean of the bounds for such periodic chains, a bound on any chain can be derived from the bounds 
on the periodic chains. 

Let us now consider general Hamiltonians possessing m types of interactions. We define a type of interaction to 
be a subset of the the local operators making up the Hamiltonian such that all the operators commute with each 
other and such that they appear with the same coupling constant in the Hamiltonian. Intuitively, we think of all 
the operators of a single type as performing the same action at different locations on the graph. For example, in the 
quantum Ising model, the set of local operators {af implementing the transverse magnetic field is considered to 
be one type of interaction. We will derive the Lieb-Robinson bound for such a Hamiltonian and provide the specific 
dependence of the Lieb-Robinson speed on the strength of the interactions. For simplicity, we will restrict ourselves 
to the case of bounded operators on a lattice. Generalisations to the commutator-bounded case and to homogeneous 
and isotropic graphs are straightforward applications of the methods shown in the previous sections, the results for 
which we state at the end of this section. 
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We consider a Hamiltonian of the following type: 

m 

H =HH h ^h ( 43 ) 

1=1 keSi 

where / labels the interaction type. We assume that all the terms in the sum are non-zero and that for alH e {1, . . . , m} 
and € Sf we have = and = ||$^|| = 1. Note that the normalization of the operators can 

be imposed without any loss of generality in the case of bounded interactions. Further, we assume that for all 
(h,h) € {1, . . . , m} 2 and i £ Si u there are exactly n^^^ j's such that j € Si 2 and T(li,i) n T(l 2 ,j) ^ 0, where n tl ^i 2 
depends only on l\ and l 2 . In addition, the graph distance travelled along the graph by passing from T{l\,i) to one 
such T(^2, j) is D^ - ^ 2 , which also depends only on l\ and l 2 . The last two properties arise from the fact that we are 
supposing the graph to be a lattice. 



A. Lieb-Robinson Bound 

For Hamiltonians with two interaction types it was easy to calculate the Lieb-Robinson bound. This is because in 
order for a signal to propagate the two types of interactions had to alternate in the operator chains which showed up 
in the bound. Therefore, we knew that half the operators were of one type of interaction while the other half were 
of the other type. For rn types of interactions, different chains could have many different combinations of each type 
of operator, the only requirement being that each operator in the chain implements a different interaction than the 
previous operator. The intuition is then to consider all possible chains which are periodic in the interaction type and 
then bound all chains by the periodic sequence of operators which has the maximal norm. This will work because 
even if a chain is not periodic, it can be subdivided into components which are subchains of periodic chains. 



Everything up to Eq.pO) does not depend on the number of interaction types. In the case of to interactions, Eq.l 20 1 
becomes 

\\K\ j (t)\\<Mj2-^r{ E 2/l ^)--- E 2Hin+i)5* n ), (44) 

where h(ik) £ {hi, h^, ■ ■ ■ , h m } and of course h(ik) 7^ h(ik+i)- This situation is now more complicated than for two 
interactions because there is not necessarily the same number of each coupling constant in each term. The power of 
a particular coupling constant in a particular term depends on the operator chain and can range anywhere from to 

\n/2\. 



In Eq.|44i the end points of the chain are separated by, at most, a graph distance of Y^k=i h 



a(ik) 



where a(i k ) e {1, 2, ... ,m} labels the operator type of i k . Hence, if d(P, Q) > ££=1 D^t* a(ih+1 \ the chain 
cannot link P to Q. As we did previously, we can bound ||-K"„V(*)II by multiplying every term in Eq. j44}by 
exp (x(J2kZl D lWT a(ik+l) ~ d ( P ' Q)))' which takes care of the end-point condition (i.e., the 5f factor in Eq.l 44 1). 
By doing this, it can be shown (see Appendix A for details) that the Lieb-Robinson bound takes the form 

\\[Op(t),O Q (0)]\\ <M(A)exp(A( 2(w ~ A 1)L(A) | t \ -d(P, Q))) , (45) 
where A is again an arbitrary positive number, M < 00 for finite A and 

L(X) = max max \ k ai ...a r e x(ai ar \l < p < q < r, (a p , a p+1 ) ^ (a q ,a q+1 ) > , (46) 

2<r<m(m— 1) (ai ,02 ,■ ■ ■ ,fl r )€{ 1,- ■ ( J 

where 



h = 

nj a\az...a q — 



» J~J h ai n ai -y ai+l , 
\ 1=1 



Eq T~i ai ~ >a i + l 
_ '=1 Uai (A7) 
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and where we assume a r+ i = a±. This implies a Lieb-Robinson speed of 

2(m- 1)£(A) 



vlr 



inf 

A>0 



A 



(48) 



We can generalise this result to homogeneous and isotropic graphs using the same methods used in the previous 
sections. All that is needed is to redefine the k's and £'s in Eq.{47}. By the definition of a homogeneous and isotropic 
graph, the result of taking the arithmetic mean of the D's ancTthe geometric mean of the n's over a chain will not 
depend on the chain taken. In other words, if we let C (P, x, 1%, Z2, ■ ■ ■ , lr) be the set of operator chains of x operators 
whose first operator's support intersects P and is of type l\ and whose q th operator is of type l g where g = [(q — 1) 
mod r] + 1, we may define: 



|C (P,a;,I 1 ,!2,...,Ir)l 



\ 



n n- 



D 



\C {P, x, Zi, Z2, 



. . , Z r ) 



EV^ jjifr+k^ijr+k + l 



(49) 



iec (p,x,h,i 2 ,.~ ,i r ) 



We now summarise our results for the Lieb-Robinson bound for a system of m interactions on a lattice or a ho- 
mogeneous and isotropic (h&i) graph. The derivation was for bounded operators but we also state the result for 
commutator bounded systems which are discussed in [17]. To state all results concisely we define a constant C and 
state the forms of fe; 1 ..j r and ^ 1 ..j r for the different possible systems: 



C 



2(m — 1) bounded, 

\f&(m — 1) commutator-bounded, 



n.i=ihkn h ^ h+1 
• VllLi h h n h i 2 ...i r 



bounded, lattice, 
bounded, h&i graph, 
commutator-bounded, lattice, 
commutator-bounded, h&i graph, 



ELi D 



D 



h...i r 



lattice, 
h&i graph. 



(50) 



(51) 



(52) 



For commutator bounded systems we define A', 



for (Zi, € {1, . . . , rn} 2 . With these redefini- 



tions, the general version of the Lieb-Robinson bound and speed are 

CL{\) 



\\[0 P (t),0, 
giving that 



v L R = inf 

A>0 



< M(A) exp ^A(- 



t\-d(P,Q)) 



CL(X) 



(53) 



(54) 



B. Calculating the Lieb-Robinson Speed 

Using Eqs. 151 1, {52} , and {54} we will calculate the Lieb-Robinson speed for a system of to interactions. If there 
are only two interactions labeled as and 1, then the max is over a set containing only one element so we have 

v LR = C inf (^ e «» lA ) = Cefeoi£oi, (55) 

A>0 [A J 

which agrees with our previous result for systems with two interactions. For more than two interactions the max 
is taken over a set containing more than one element which brings about a qualitative difference: One obtains a 
different functional dependence on the strengths of the interactions (the h's), as well as the structure of the graph 
(the n's and D's) as compared to the algebraic dependence in the case of two interactions. 
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To calculate vlr for more than two interactions, we must take the maximum of a set of functions of the form 
fe»i...< r e 1, " ir • Since the maximal element of this set can change for different values of A, some care must be taken. 
In view of Eq.|54| we divide these elements by A since it does not change which element is maximal and we call this 
set A. That is 

■■ {•"■•-•;• (56) 



for the different possible subchain labellings i± . . . i r , as given in Eq.|46 1. Further, we define a function A mllx (X) — 
max A. Therefore 

v LR = Cmf o A max (\). (57) 

The element of A which has the largest value of k will be maximal for A — > 0. The value of A at which two elements 
of A intersect is given by 1 - r A = J1 " " Js A and corresponds to the value 

xCH...i*> _ M fc 'i...* t ) /cov 
A (fi-i.) 7 T7 ■ 

sji-.-is Mi- ..it 

Notice that for A^'''^, to be positive, we must have £, n ...j s > £ii...z«/ which is also the condition for the two curves 
to intersect. Now suppose that k^...^ is the largest value of k of the elements in A and that A^"'|\ is the smallest 
intersection point involving this element. Next, suppose that A[^ ^ is the smallest intersection point greater than 
involving the element labeled by ji . . . j s . Continuing in this way, the function A max (X) will be given by 

< A < A^- 3s) 



(A) = I fcj -"V" > A|i-;tJ<A<A^, (59) 







.i A 


kjl ■ 


A 






A 




™J1 ■ 







r) ' 



Since the function — " 3a A has only one minimum at A = l/^j 1 ...j s , it follows that if X\" '^j < 1/Cji...j s < 



Ayj.',.^) for any one of the intervals, then the infimum of A max (A) is given by -4(1/^...^). Otherwise, the infimum of 



AmnxW will be at one of the points A^ '/^ . This completes the algorithm for finding the Lieb-Robinson speed v LR 
for any finite number of interactions. 

For concreteness we work out some general scenarios for the case of three interactions. We label the interactions 
by 1,2, and 3 in which case 

ffci 2 e« 12A fc 23 e« 23A k 31 eA 3lX k 123 e^ 123X k 132 e^ 132X 



AAA A A 

,£123132-*- I 

(60) 



fc 213ie 6i3iA fci 232 e^ 32A fc 1323 e^ 23A fci 2 3132^ 123132A 

A ' A ' A A 



Note that we have only considered distinct elements of the set. For example, we did not consider fc23ie A 231 because 

A A 

There are two possible scenarios. The first is that one of the k's and its associated £ are greater than all the other 
k's and £'s. In that case, one element of A will be greater than all the others regardless of the value of A and the 
Lieb-Robinson speed will be given by minimising that particular element of A. If, for example £132 is greater or 

equal to all the other fc's while simultaneously £132 is greater or equal to all the other £'s, then fcl32e A 132 will be the 
maximum element irrespective of A, and so we have 

n . t fci 32 e 5l32A 

vlr = C mf , 

A>0 A 

= Cek 132 £, 132 . (61) 
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TABLE I: The minimal point of the function A max (A) is located at A = A^f . The top row indicates whether the function j4 max (A) 
is increasing or decreasing in the intervals of A determined by the various critical points. The middle row indicates the maximal 
element of the set A while the bottom row is a scale for A beginning at A = 0. Also shown are f^ 1 , £23 / an d $123 which are the 



minima of their respective elements in A. In this case they are not minima of A n 
which their corresponding element of A is maximal. 



( (A) because they do not lie in the intervals for 



For a bounded system on a lattice C = 2(m — 1) and 



^123 = \Jhih 2 h 3 ^/ni^ 3 n 3 ^ 2 n 2 ^,\, 

_ D 1 ^ 3 + D 3 ^ 2 + D 2 ^ 1 
£123 = 7, 1 



(62) 
(63) 



which allows us to solve for the Lieb-Robinson speed explicitly in terms of the coupling constants and lattice dis- 
tances as 



4e 



V L R = y \/^^M^T(0^ 3 + D 3 ^ 2 + D 2 ^ 1 ). 

In the second scenario, no single element of A is greater than all the others for all values of 
that case, that we have k 12 > k 23 > k 123 > k 31 > k l32 > k 1232 > £123132 > ^2131 > ^1323 an 
£1323 > 6131 > 623132 > £1232 > £31 > £i32- With these assumptions one finds that 



(64) 

of A. Let us suppose, in 

and 623 > £23 > €12 > 



c(A) = 



( k 12 e 6 ^ x 
X ■ 

X 



fcl23C 



X 



0<A<AS, 



\ <23) 

A (12) 
\ (123) 
A (23) 



< \ < X <123) 

<A<..., 



(65) 



and also that £^ , £^ , an d £123 do not lie in the intervals of A for which their corresponding element in A is maximal, 
i.e. none of these are the minimal point of j4 max (A). 

In Table U we state the maximal element of A in the different relevant intervals of A. We also indicate whether 
the function j4 max (A) is increasing or decreasing in each interval. From the table we see that the function ^4 max (A) is 
minimized for A = X[Hf . 

Thus we have 



VLR = C 



k 23 e 



5, (123) 
A (23) 



(66) 



For a bounded system on a lattice we have C = 2(m — 1) and 



^23 = V 'h 2 h 3 ^/n 2 ^ 3 n 3 ^ 2 , 
D 2 ^ 3 + D 3 ^ 2 



523 



(67) 
(68) 
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which allows us to solve for the Lieb-Robinson speed explicitly in terms of the coupling constants and lattice dis- 
tances as 

v LR = 2(m - 1) ™> r fc| 3 123 - fa3 fc«|r ei23 (69) 

[ h 2 + h 3 _ h 1 +h 2 +fl3 1 , r w 2 ->3+"3->2 _ »1^2+»2^3+"3-»l 1 

_ A 2 3 LZJ 2 3 1 

J>l->2 + J>2->3 + J>3-H _ g 2->3 + rj3->2 

3 2 

g 1 ~ f2 ^ z? 2-»3^. rJ 3— »1 



/- £jl -^-2 + D 2^.:i j7_> — > l r>2— >3 + j>j^2 

x V n 2 n 3 n2^3n 3 ^2 3 2 

^2— >3 + z -)3— >2 

2 

3/, , , gggg+pggj l->2 + D 2^3 +r) 3-H 

x V h 1 h 2 h 3 n 1 ^ 2 n2^n- i ^ 1 2 3 . (70) 

In general, when the minimum of A max is attained at a minimum of one of the elements in A, the Lieb-Robinson 
speed will be of the form of Eq.[|6"l"). When the minimum of A max occurs at the point where two elements of A are 
equal, the Lieb-Robinson speed will be of the form of Eq.([69). Thus, the Lieb-Robinson speed will always be of one 
of the two forms 

v s LR (k,0 = 2ce(m-m (71) 

or 

12 ei 

■k(h.h,6.6) = -fe) - m£T • <72) 

where c = 1 for bounded systems and c = \/2 for commutator-bounded systems. This is true also for any m > 3 since 
adding more interactions merely increases the number of elements of A. Now that we have seen how to calculate 
the Lieb-Robinson speed for any number of interactions, let us consider a physical example. 

C. Example: XY-Model 

We investigate the anisotropic quantum XY-model on a 2-dimensional square lattice. This model is well under- 
stood and is studied mainly in connection with quantum phase transitions. This model is also of practical interest, 
for example, in understanding the isolator-superconductor phase transition in cuprate type II superconductors 1 23 1 . 
At zero temperature the system exhibits two phases, one with an energy gap and another which is gapless. 

The Hamiltonian for the quantum XY-model is given by 

H = J ( S n S ?n + S y n Syj + Dj2(S z n f, (73) 

(n,m) n 

where (n, m) represents the sum over nearest neighbours on the sites, n, of a square lattice. Note that the gapless 
and gapped phases are determined by S < S c and 5 > S c respectively where S = and where 5 C has been estimated 
to be approximately 3.6 [24 1. For S > 5 C , the correlation length of the system approaches a finite value leading to a 
quantum paramagnetic ground state with no long-range order. We label the S X S X , S y S y , and S z operator types by 
X, Y, and Z respectively 2 . 

Since the fc's and £'s are symmetric under cyclic permutations as well as the interchange of X and Y, the set A will 
have only four distinct elements: 

f k XY e x ^ k xz e x ^z kxyze KxYz kxYZY e x ^v^ \ 

1 A ' A A ' A I' ( ' 



2 Note that even though S% and S„ have the same coupling constant, we must nevertheless consider them as two separate interactions 
because otherwise the conditions laid out between Eq.|ll| and Eq.|12| would not be satisfied. 
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Note that XZYZ = XZ and XYZXZY = XYZ by breaking the larger chains in two and using the symmetries. 

The operators X and Y are located on edges of the graph, whereas the operators of type Z are located on the 
plaquettes, the elementary squares of the lattice composed of four edges. Moving from a plaquette to an edge 
corresponds to moving a lattice distance of 1 /2, while moving from a plaquette to another plaquette corresponds to 
moving a lattice distance of 1 . Therefore 

D x ^z = D Y ^ Z = 1/2, (75) 
Dx^y = 1. (76) 

Note that these distances are symmetric. Furthermore, from an edge one can move to one of two plaquettes or one 
of 7 other edges (including the one currently on). From a plaquette one can move to one of four edges. Therefore, 

n x ^z = n Y ^z = 2, (77) 
n Z ->x = n z ^Y = 4, (78) 
n x ^Y = n Y ^x = 7. (79) 

Using this information and the coupling constants hx — J, hy = J, and hz = D in Eqs. ( |5"T) and | |52~) , we get 

k X y = 7JS, fxr = 1, 

k xz = V&JDS, txz = \, 

kxYZ = {/(7 ■ 8)J 2 DS, txYZ = | 

kxYZY = </(7 2 • 8) J 3 DS, Uyzy = \, (80) 

where S is a half integer giving us the type of spins on the lattice. From Eq.j58), we obtain 

\{XY) = \(XY) = \(XY) — \{XYZ) — \(XYZY) _ \(XYZY) — ^ lll(4/7) + 111 f = ] . (81) 

(XZ) {XYZ) (XYZY) {XZ) {XZ) {XYZ) \ I / • I ^ j I \ I 



This implies that all of the elements of A intersect at the same value of A. First note that this value is negative 

A 



for D/2J < (7/4) 2 and thus the curves do not intersect. In this case, kxYe XY is maximal for all A and so the 



Lieb-Robinson speed is 

v LR = 2e(m - l)k X yCxY, 

= 28eJS. (82) 
For D/2J > (7/4) 2 we have 



k X Ye SxY X 



< A < A< 



Afsrv, < A < 00. 



A ' '\XY) 



First let us consider the cases when the minima occur within the intervals. The minimum occurs in the first interval 
when < £ XY < \\xy\, which corresponds to (7/4) 2 < D/2J < e(7/4) 2 and gives the same speed as for D/2J < 
(7/4) 2 . The second interval contains the minimum when \{ X y) < £xz ^ °°/ which occurs for D/2J > e 2 (7/4) 2 . In 
this case the minimum value is A max (£ xz ) and gives a Lieb-Robinson speed of 

v LR = 2e(m - \)k X zixz 

= 4eV2JDS. (84) 
Finally, in the region e(7/4) 2 < D/2J < e 2 (7/4) 2 the minimum occurs for A = \\ XY \ and so 
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V L R 


Conditions 


Simplified Cond. 


v S LR(kxz,S,xz) 
= 4eV2JDS 


x< Xi, > > f- 1 

"(JZ) — >>xz 

\iXYZ) v, ^-1 
A (XZ) — SXZ 

,(xyzv) ^ ^-1 

"(JZ) — SJ!fZ 


&>(J)V 


VLR(kxz,k X Y, 
£xz,£xy) 

32 DS 
- 7 21n(f) + ln( ft) 


j 1 , ( (VV| 

sxy — A (xz> 

x< Xi, > < t- 1 

A (xz) — Sxz 

A (XYZY) — A (XZ) 

A (xyz) — A (xz) 


(I) e<£<(|) e 2 


V S LR (k X Y,ZxY) 

= 2SeJS 


x< xv ') <r f- 1 
A (xrz) — sxy 

x< Xi ') < t-l 

A (IZ) — S1Y 
\(") < e-1 





TABLE H: The Lieb-Robinson speed for the XY-model. The first column gives the functional form of the Lieb-Robinson speed in 
terms of the interaction strengths. The second column gives the region of applicability of each of the three functional forms of the 
Lieb-Robinson speed in terms of the A's and £'s. The last column gives the regions in terms conditions of the interaction strengths. 
The results of this table are depicted graphically in Fig. [T] 



T(XZ) 
A (XY) 



7 21n(4/7) + ln(Z>/2J)' v "" ; 

The results and the various conditions are summarized in Table {P) . The first column gives the Lieb-Robinson 
speed in the region of phase space defined by the inequalities of the third column. By phase space, we simply mean 
the space M. 2 ^ of all possible values of the coupling constants J and D. The different values of the Lieb-Robinson 
speed in the different regions of the phase space are plotted in Fig. [T] In general, the number of functional forms 
for the Lieb-Robinson speed could be as many as twice the number ordistinct intersection points obtained in Eq. (jHTJl 
plus one, i.e. the minimum could occur at each intersection point or in each interval. 

Also plotted in Fig. [T|is the critical line for the quantum phase transition (dotted line) which occurs for sw 3.6. 
It is important to remember that the regions of phase space with different Lieb-Robinson speeds do not correspond 
to different quantum phases. Quantum phases occur in the ground state at a particular point of the phase space, 
whereas the Lieb-Robinson speed is the maximum speed of any signal in any state and not just small excitations 
over the ground state. 



v LR = 2(m - 
32 



Conclusions 



In this article we showed how the Lieb-Robinson bound could be considerably tightened for the case of a finite 
number of interactions. We also derived a bound which was sensitive to the relative interaction strengths. We 
first considered the case of two interactions, which is one of the most frequently encountered cases. Indeed, for a 
physically relevant bound we must use not the full Hamiltonian, but the effective Hamiltonian near the ground state, 
which often contains only two interactions. We then generalised the result to any finite number of interactions. We 
observed that for m > 3 interactions, the space of coupling constants was subdivided into a finite number of regions, 
each of which produced a different functional form for the Lieb-Robinson speed. Interestingly, the dependence 
of the speed on the coupling is not necessarily algebraic anymore for three or more interactions as it was for two 
interactions; it was found that the logarithm of the coupling constants could appear in the speed. 

It is also important, especially for applications in such domains as quantum gravity, that we be able to calculate 
the bound not only for a system on a lattice but also for general graphs. We derived explicit equations for the Lieb- 
Robinson bound and speed in the case of homogeneous and isotropic graphs. This is the best we can hope to have 
since if the graph is not homogeneous and isotropic, the speed will not be constant. 
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FIG. 1: The (DJ)-plane of the XY-Model and the three different formulas for the Lieb-Robinson speed. If > (|) 2 e 2 / then the 
Lieb-Robinson speed is 4e\/2JDS; if (|) 2 e < ^ < (|) 2 e 2 , then the Lieb-Robinson speed is T 2 i n ( 4 D ) ' anc ^ nna %' ^ 

^7 < (i) 2 e ' tne Lieb-Robinson speed is 28eJS. The dotted line indicates the critical line for the phase transition which runs 
along the line ^7 = 3.6 



Appendix A: Multiple Interactions Calculations 

l—l n a(ifc)— >a(i fc + i) 



In Eq.|44 1 the end points of the chain are separated by a graph distance of J2k=i ^ a (il) k+1> . Hence, if d(P, Q) > 

Z)fc=i ^f(it)^ ^ fe+1 ^' tne chain cannot link P to Q, where a(i k ) G {1,2,..., to} labels the operator type of i k and 
/i(ifc) € , hz,..., h m } gives the coupling constant of operator ik- We can thus legitimately multiply every term by 

cxp (x( J2t=l D a a tn!r a{ik+l] ~ d ( P > < 3))) and dro P the end -p° int condition. This gives us: 



<_XJ I i n ft— J- 

|^VWII<M^^L( J2 2h(i 2 )... Y: 2h( ln+1 ))ex P (\(Y,D:^ a ^-d(PQ)) 

n=Q ' i2&Z il i n+1 eZ in k=l 

^M^^±^{Y.^ D ^ r "--- E M^ +1 )e<-^ +l) ). (86) 



n=0 i2eZ il i n+1 £Zi, 
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The equivalent equation for commutator bounded systems is: 

\\<\ j (t)\\ 



- •t 



^ M E^r( E 2 %)%).w- E 2^ a(in)a(i „ +0 / l (i„ + o)e X p(A(^^;^ +l) -d(p ) g)) 

n=0 ' j 2 e2 tl in+i£Z(„ fe=i 

(87) 

For (I1J2) G {1, . . . , m} 2 , ii 7^ Z 2 and A > 0, we define 



Ni t i 2 (A) = ni 1 ^i 2 hi 2 e XDl i for bounded systems 
iVj^j (A) = ni 1 ^i 2 Ki 1 i 2 hi 2 e XDl i for commutator-bounded systems 

N hll (X)=0, (88) 

We proceed for bounded systems but note that the commutator-bounded derivation is virtually identical. Now, 
from an operator of type 1, there are rix-yi ways to choose an operator of type I as the next in the chain and there are 
771—1 possibilities for I € {2, . . . , m}, and similarly for the other operator types. So, if we consider the embedded 
sums 

E Hh)e XD «n) ... £ ^(*n+i)e AD »<-) , (89) 

we see that each individual sum of the type 

E ft (^+i) e (90) 
can be written as 

E <!S"V AD -^ ! = e ^mW. (9i) 

This means for fixed n we can write 

EM* 2 )e A <^ a<l2, ...E Kin + i)e A ^' H ' ( '" +ll = E ^a(n)a 2 (A)... E ^„a„ +1 (A). (92) 



of Eq.|86 , where the ai's represent interaction types. We thus have that Eq.|92 1 has (m — 1)" terms, each of which 
has theform 

N ai a 2 (X)N a2a3 (X)N a3ai (X) . ..N anan+1 (X). (93) 

Let us call a linking sequence ofN's of length q a finite sequence of 2Vy (A)'s TV ai6l (A) , N a2ll2 (A) , . . . , N aqbq+1 (A) such that 
a; + i = b[. Considering that there are only to interactions, we have that the set of transition factors {JV a (,(A)|l < a < 
to, 1 < b < m,a ^ b} contains to(to — 1) elements; thus, any sequence of m(m — 1) + 1 TV's will necessarily have at 
least one repetition. Hence, there are a finite number of linking sequences of TV's such that all TV's in the sequence are 
different, and they are of length q < m(m — 1). Let us call a linking sequence of TV's such that all TV's in the sequence 
are different and such that a\ = b q+ \ an elementary cycle. To each elementary cycle TV 0ia2 (A), TV a2a3 (A), . . . , TV a?ai (A) 
we associate two numbers, 



h 



\ 1=1 



- '=i Uai (94) 
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where we have made the identification ai = a g +i. We then have that 

iV aia2 (A)7V a2a3 (A)...iV Q?ai (A) = (fca lBa ...«,e^«--«)*- (95) 



Now Eq.|93 i is the product of a linking sequence of N's of length n. In the first m(m — 1) + 1 ./V's of the product there 
must be, as previously mentioned, at least two identical TV's; let us now call them N bl b 2 . Because we are dealing with 
the product of a linking sequence of N's, the subsequence of N's starting with the first N blb2 and ending with N just 
before the second N blb . 2 is an elementary cycle of the type N blb2 , N b2ba , . . . , N bqbl and if we remove it from the full 



sequence, the remaining sequence is still a linking sequence of length n — q. So we can now write Eq.(93 1 as 



(fc 6l62 ... fc ,e A ^ 6 -^) 9 iV ClC2 (A)...iV c „_ 9C „ +1 _ 9 (A), (96) 

where iV ClC2 , . . . , N Cn _ qCn+1 _ q (A) is a linking sequence of N's. This process can be reiterated until at most m(m — 1) 
TV's are left. 



^13 ^32^24^43 ^32 ^24^42 ^24^43^32^21^13^34 

(fc 3 24e^324)3 (fc 24 e*«24)2 

= (W AC324 ) 3 (fc 24 e A ^) 2 jVi3 ^32^24^43 ^32^21^13^34 

(fc 324 e^S324)3 

= (fc324e Afa4 ) 6 (fc24e AC24 ) 2 jVl 3 iV32jV2l iVl3^34 

(fc 132 e^i32)3 

= (fc324e A « 324 ) 6 (fc24e A « 24 ) 2 (fci32e A « 132 ) 3 ^Vi 3 iV34. 

(97) 

Eq. |97)is an example of the rewriting of a term of 
the form in Eq.|93) as a product of the elementary 
cycles defined m~Eq.l|94). The dependence on A 
is not written explicitly in the above equation for 
simplicity. 

So if we now define 

L(X)= max max < k ai ... ar e xia ^ ar \l < p < q < r,(a p ,a p+1 ) ^ (a q ,a q+1 )[ (98) 

2<r<m(m— l)(ai,a2,... : a r )€{l,...,m} r ^ J 

G(A)= max , I ^^-\N lzU {\) ± o) < oo, (99) 



(Ji,WsA)e{i,...,m} 4 {N hu (X) 
we conclude from the previous argument that 

iV aia2 (A)7V a2a3 (A)7V a3Q4 (A)...7V a „ Q „ +1 (A) < GiXy^-'Hmr. (100) 



Since there are (m — 1)™ such terms in Eq.|92 1, this in turn gives us that 



J2 h(i 2 )e X <^ a(t2) Ki n+1 )e XD t^ a{ ' n+l) < G(X) mt - m ~^ ((m - l)L(A))" . (101) 



Inserting this last inequality in Eq.|[86), we obtain 

°° I Of I™ 

= M(A) cxp (A( 2(m ~ A 1)£(A) | * | -d(P, Q))) , (102) 
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where M (A) = M(G(A)) M ( m_1 ). This last inequality implies a Lieb-Robinson speed of 2(m * )£(A) for any positive A. 
To obtain the tightest bound, we take the infimum over A 

. r 2(m- l)L(A) Mnn . 

vlr = inf ^ r^^- 103 

A>0 A 
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